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We present the supersymmetric Witten and double Darboux (strictly isospectral) constructions as applied to the diffusion 
of thermal neutrons from an infinitely long line source. While the Witten construction is just a mathematical scheme, the 
double Darboux method introduces a one-parameter family of diffusion solutions which are strictly isospectral to the stationary 
solution. They correspond to a Darboux-transformed diffusion length which is flux-dependent. 
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Supersymmetric one-dimensional (ID) quantum me- 
chanics has been considered by Witten in 1981 as a toy 
model for symmetry-breaking phenomena in quantum 
field theory [Q . With great speed its status has changed 
to a powerful research area as one can contemplate in the 
most recent review [|[. In a series of papers, we obtained 
interesting results for various physical problems by us- 
ing Witten's factorization procedure of the Schrodinger 
ID operator and a more general supersymmetric double 
Darboux method |^ introduced by Mielnik [Q. The aim 
of this work is to apply these two simple, supersymmetric 
methods to the theory of diffusion of thermal neutrons. 
We shall use the illustrative example of the neutron dif- 
fusion problem as presented in the textbook of Arfken ||] 
referring to an infinite line (Dirac delta) source of neu- 
trons, thus providing the Green's function for this case. 
The steady state continuity equation for the neutrons 
reads 



^ + 5 = 



(1) 



where the first term represents the diffusion, the second 
stands for the absorption losses and S stands for the 
source strength. The diffusion constant D is related to 
the neutron mean free path as follows D — ■j^^i_^2/3A) ' 
where A is the atomic number of the scattering nucleus 
and enters as a correction for the anisotropy of scatter- 
ing in the laboratory system. The solution (p has the 
physical meaning of neutron flux being the product of 
neutron density times average velocity. Finally, Sa is 
the macroscopic absorption cross section, i.e., the prod- 
uct of the microscopic (atomic) absorption cross section 
and the number of atoms per unit volume. Usually it 
is assumed that the absorption is small compared to the 
scattering. As we said, the neutron source is considered 
as an infinitely long line source that can be taken along 
the z-axis embedded in an infinite diffusing medium. Its 
strength is neutrons per unit length per unit time. 
Thus S — So6{p) where S{p) is the cylindrical Dirac delta 
function. For p ^ 0, in view of the no z and 9 depen- 
dence, one gets the radial diffusion equation 



dp'^ ~^ ^ dp 



(2) 



where we have made use of the "diffusion length" k^^ — 
a/ D /Yia- Then one can write the general solution in 
terms of the modified Bessel functions </> — aiIo{kp) + 
a2Ko{kp). The physical solution is only the Kq term be- 
cause the flux is supposed to decrease at large distances. 
The constant 02 may be determined by requiring that D 
times the integral of the negative gradient of the neutron 
flux around the edge of a small pillbox of unit height be 
equal to the production strength Sq within the pillbox. 
The box is small (p — )■ 0) to eliminate absorption and 
gives 



Sq — lim Da2 



'WKoikp) ■ po]pdd , 



which is a two-dimensional form of Gauss's law. 
the series form of Ko{kp) this turns into 



So 
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Using 



(4) 



(5) 



for the final form of the common solution. 

Let us pass now to the supersymmetric constructions. 
We have already presented these schemes in previous pa- 
pers H. One needs a self-adjoint form of Eq. (2) that 



can be obtained by the change of function 
leading to 



-1/2 



^ -(fc2 



(6) 



If Eq. (6) is interpreted as a Schrodinger equation at fixed 
zero energy with the potential Vb[p) = fc^ — then 
one can think of supersymmetric quantum mechanical 
methods as follows. 
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(i) . Witten's construction has two steps. The first 
is the factorization of Eq. (6) by means of the opera- 
tors Ai = + W{p) and ^2 = ^ - W{p), where W 
is the so-ca ifed superpotential function that can be de- 
termined from the initial ( "bosonic" ) Riccati equation 
Vb — — , or even more directly as the negative 
of the logarithmic derivative of the stationary solution ip 
{W — — ^In-!/;). The second step of the Witten con- 
struction means to pass to a "fermionic" problem by 
merely changing the sign of the derivative term in the 
Riccati equation and determining the new potential Vp, 
which is interpreted as a "fermionic" partner of the initial 
potential. Thus, Vp = + W . This partner potential 
enters a "fermionic" equation for which the factoring op- 
erators are applied in reversed order. On the other hand, 
the superpotential should be the negative of the loga- 
rithmic derivative of a "nodeless" solution, which is our 
case since modified Bessel functions of zero order occur. 
Moreover, the singularity at the origin of the superpoten- 
tial function is not disturbing because we have an infinite 
line source there. 

(ii) . The double Darboux construction in the super- 
symmetric framework allows one to use the general so- 
lution Wgen of the "bosonic" Riccati equation and not 
the particular one as for the Witten construction. Thus, 
the factorization operators are now Ai = ^ + Wgen and 



Wa. 



In this way one can introduce a one- 



_d 

dp ' ' gen 

parameter family of "bosonic" potentials having the same 
"fermionic" partner, i.e. 



— \n{I{p) + X) 



which is a general Darboux transform of Vb ■ Eq. (7) can 
be written in the following form 



vUp,tA) = VB{p)- 



I{p) + A + {I{p) + Xf 



(8) 



where I{p) = ip'^{r)dr and A is the family parameter, 
which is a real positive quantity, being as a matter of 
fact the Riccati integration constant ||^. Besides, there 
is a modulational damping of the general solution, which 
reads 



0(p; A)- v/A(A + l) 



lip) + A ' 



(9) 



where the square root is a normalization constant. One 



can check easily that Wgen = 



-iinVXp;A). 



Some 



plots of the strictly isospectral potentials Eq. (8) and 
the isospectral solutions Eq. (9) are presented in Figures 
1 and 2 for different A values and Oi =02 = 1; the 
case A = 00 corresponds to the original "bosonic" solu- 
tions. The interesting point is that one can work with 
the general flux solution and still get physically mean- 
ingful, locahzed-type solutions (see Fig. 3). We mention 



that in quantum mechanics irregular wavefunctions have 
been found useful in state reconstruction 0. To see the 
physical meaning of the isospectral construction for neu- 
tron diffusion we make the change of function ip = ^/pcj) 
in the Schrodinger equation ip — Visoip = to get 



P 



dp^ 



P- 



where k^^^ 



dp 

= [fc2 - 2^ \n{I{p) + A)] -1/2 is an effective, 
flux-dependent diffusion length, which is the Darboux 
transform of the constant diffusion length. Thus, one can 
see the connection between the common radial neutron 
diffusion equation and a similar isospectral diffusion. 

In conclusion, we have shown that simple supersym- 
metric methods of ID quantum mechanics may provide 
interesting results in neutron diffusion. 
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Figure captions 

Fig. 1: The original potential Vb = fc^ — ^ 
fc = 1 and four of the strictly isospectral potentials for A 
taking the values 1, 1000, 3000, 6000, respectively (the 
corresponding wells are from left to right, respectively). 

Fig. 2: The original (nonphysical) solution tp for the 
superposition constants ai = 02 = 1 and the strictly 
isospectral solutions for the same values of A as in Fig. 
1. 

Fig. 3: The flux solutions (j> corresponding to the tp 
ones in Fig. 2. Unpublished note: In an inverted poten- 
tial the isospectral 'flux' solutions could be interpreted 
as a sort of 'resonant' states. 
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Fig. 2 
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